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Abstract. Introduced in 1963, Glauber dynamics is one of the most 
practiced and extensively studied methods for sampling the Ising model 
on lattices. It is well known that at high temperatures, the time it takes 
this chain to mix in on a system of size n is O(logn). Whether in 
this regime there is cutoff, i.e. a sharp transition in the L^-convergence 
to equihbrium, is a fundamental open problem: If so, as conjectured 
by Peres, it would imply that mixing occurs abruptly at (c + o(l)) logn 
for some fixed c > 0, thus providing a rigorous stopping rule for this 
MCMC sampler. However, obtaining the precise asymptotics of the 
mixing and proving cutoff can be extremely challenging even for fairly 
simple Markov chains. Already for the one-dimensional Ising model, 
showing cutoff is a longstanding open problem. 

We settle the above by establishing cutoff and its location at the 
high temperature regime of the Ising model on the lattice with periodic 
boundary conditions. Our results hold for any dimension and at any 
temperature where there is strong spatial mixing: For this carries all 
the way to the critical temperature. Specifically, for fixed d > 1, the 
continuous-time Glauber dynamics for the Ising model on (Z/nZ)'' with 
periodic boundary conditions has cutoff at (ci/2Aoo) logn, where Aoo is 
the spectral gap of the dynamics on the infinite-volume lattice. To our 
knowledge, this is the first time whore cutoff is shown for a Markov chain 
where even understanding its stationary distribution is limited. 

The proof hinges on a new technique for translating L^-mixing to 
//■^-mixing of projections of the chain, which enables the application of 
logarithmic-Sobolev inequalities. The technique is general and carries 
to other monotone and anti-monotone spin-systems, e.g. gas hard-core, 
Potts, anti-ferromagentic Ising, arbitrary boundary conditions, etc. 

1. Introduction 

The total- variation cutoff phenomenon describes a sharp transition in the 
L^-niixing of a finite ergodic Markov chain: Over a negligible time period, 
the distance of the chain from equilibrium drops abruptly from near its 
maximum to near 0. Though believed to be widespread, including many 
important families of chains arising from statistical physics, cutoff has been 
rigorously shown only in relatively few cases (ones where the stationary dis- 
tribution is completely understood and has many symmetries, e.g. uniform 
on the symmetric group) . Here we establish cutoff for Glauber dynamics for 
the Ising model, one of the most studied models in mathematical physics. 
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Already establishing the order of the mixing time is in many cases chal- 
lenging, with an entire industry devoted to the study of such problems. 
Proving cutoff and its location entails not only obtaining the order, but also 
deriving the precise asymptotics of the time it takes the chain to mix. In his 
1995 survey of the cutoff phenomenon Diaconis [7] wrote "^t present writing, 
proof of a cutoff is a difficult, delicate affair, requiring detailed knowledge 
of the chain, such as all eigenvalues and eigenvectors. Most of the exam- 
ples where this can be pushed through arise from random walk on groups, 
with the walk having a fair amount of symmetry" . To this date this essen- 
tially remained to be the case, and present technology (e.g., representation 
theory, spectral theory, techniques for analyzing 1-dimensional chains, etc.) 
does not suffice for proving cutoff in high- dimensional chains with limited 
understanding of their stationary distribution, such as Glauber dynamics 
for the Ising model (stochastic Ising model) on the 3-dimensional lattice. 

Introduced in 1963 |16 , Glauber dynamics for the Ising model on the 
lattice (see Section |1.1| for formal definitions) is one of the most practiced 
methods to sample the Gibbs distribution, and an extensively studied dy- 
namical system in itself, having a rich interplay of properties with the static 
stationary distribution. For instance, as we describe in Section |1.1[ it is 
known that on (Z/Z„)^, at the critical inverse-temperature /3c for uniqueness 
of the static Gibbs distribution, the spectral gap of the Markov semigroup 
generator of the Glauber dynamics exhibits a phase-transition from being 
uniformly bounded to tending to exponentially fast in n. It is further 
known that at high temperatures on (Z/Z„)'^ this dynamics mixes in time 
O(logn), yet the precise asymptotics of the L^-mixing time were unknown 
even in the one-dimensional case: It is an open problem of Peres (cf. [2 3 J 
to determine cutoff for the Ising model on 

The only underlying geometry for which cutoff for the Ising model had so 
far been established is the complete graph (|13,22 ), where the high symme- 
try reduces the analysis to a birth-and-death magnetization chain. However, 
this sheds no light on the existence of cutoff for lattices, where there is no 
such reduction. Peres ( p2|[23] ) conjectured that in any dimension d, Glauber 
dynamics for the Ising model on (Z/Z„)'^ should exhibit cutoff. 

Our main results confirm the above conjecture and moreover establish 
cutoff and its location in a wide range of spin system models and geometries. 
We first formulate this for the classical two-dimensional Ising model. 
Theorem 1. Let (3c = ^log(l -|- \/2) be the critical inverse-temperature for 
the Ising model on I?. Then the continuous-time Glauber dynamics for the 
Ising model on {'L/n'L)'^ at inverse-temperature < P < Pc with periodic 
boundary conditions has cutoff at X^logn with a window o/ O(loglogn), 
where Aqo is the spectral gap of the dynamics on the infinite-volume lattice. 
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Figure 1. Cutoff plienomenon for the (total- variation) 
distance from stationarity along time in Glauber dynamics 
for the Ising model on Z^, as established by Theorem [2| 
Highlighted region denotes the cutoff window of O (log log n). 

In the above theorem, the term cutoff window refers to the rate at which 
the L^-distance from stationarity drops from near 1 to near 0. More pre- 
cisely, let tMix(£) be the minimum t > where the heat-kernel Ht associated 
with a Markov chain is within a total-variation distance of e from station- 
arity. A family of chains is said to exhibit cutoff if for every fixed < e < | 
we have tMix(£)/iMix(l — e) ^ 1 as the system size tends to oo. A sequence 
Wn is said to be a cutoff window if t^nxi^s) = tMix(l — ^) + 0{wn) for every e. 

Our results hold for ifLjriX)'^ in any dimension d, inverse-temperature 
/? and external field h so that the corresponding static Gibbs distribution 
has a spatial dependence property known as strong spatial mixing (and also 
as regular complete analyticity) . This property, defined by Martinelli and 
Olivieri in their seminal paper [30], holds in all regimes where O(logn) 
mixing is known for Glauber dynamics for the Ising model. In particular 



(see 30 - 32 ) , on there is strong spatial mixing for any (3 with an external 
field 7^ 0, as well as for any < /3 < /3c when h = 0. 

The next result settles the conjecture of Peres for cutoff for the high 
temperature regime of the Ising model in any dimension d> 1. 

Theorem 2. Let d> 1 and consider the continuous-time Glauber dynamics 
for the ferromagnetic Ising model on (Z/nZ)'^ with periodic boundary con- 
ditions, inverse-temperature (3 and external field h. Suppose that P, h are 
such that there is strong spatial mixing. Then the dynamics exhibits cutoff 
at {d/2Xoo) logn with a window o/0(loglogn), where Aqo is the spectral gap 
of the dynamics on the infinite-volume lattice. 
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In the special case of d = 1 and no external field, it is known that strong 
spatial mixing always holds and that the spectral gap at inverse-temperature 
P is 1 — tanh(2/3) independent of the system size (cf., e.g., [23]). We thus 
have the following corollary to establish the asymptotic mixing time of the 
one-dimensional Ising model, answering the aforementioned question of [23] . 

Corollary 3. For any /3 > 0, the continuous-time Glauber dynamics for the 
Ising model onTLjuTL with 'periodic boundary conditions, inverse-temperature 
13 and no external field has cutoff at |(1 — tanh(2/3))~-^ logn. 

Our proofs determine the cutoff location in terms of A(r), the spectral gap 
of the dynamics on the d-dimensional lattice (Z/rZ)'^ for a certain r = r(n). 
As a biproduct, we are able to relate the spectral gap on tori of varying sizes 
and obtain that they converge polynomially fast to the spectral gap of the 
dynamics on the infinite-volume lattice. 

Theorem 4. For d > 1 let A(n) be the spectral gap of the continuous-time 
Glauber dynamics for the Ising model on (Z/nZ)"^ with inverse-temperature 
P >0 and external field h. If there is strong spatial mixing for h then 

|A(n)-Aoo|<n-i/2+°W, 

where Aqo is the spectral gap of the dynamics on the infinite-volume lattice. 

1.1. Background and previous work. While our results hold in greater 
generality, we will focus on single-site uniform interactions for the sake of the 
exposition: The Ising model on a finite graph with vertex-set V and edge- 
set E is defined as follows. Its set of possible configurations is O = {±1}^, 
where each configuration corresponds to an assignment of plus / minus spins 
to the sites in V . The probability that the system is in a configuration cr G $7 
is given by the Gibbs distribution 

^(^) = ^ exp ( /3 ^ a{u)<j{v) + h^, <u) J , (1.1) 

\ uveE ueV / 

where the partition function Z(P) is a normalizing constant. The param- 
eters P and h are the inverse-temperature and external field respectively; 
for /3 > we say that the model is ferromagnetic, otherwise it is anti- 
ferromagnetic. These definitions extend to infinite locally finite graphs (see 
e.g. (24||28]). 

We denote the boundary of a set A C ^ as the neighboring sites of A 
in y \ A and call r G {±1}^^ a boundary condition. A periodic boundary 
condition on (Z/nZ)'^ corresponds to a d-dimensional torus of side-length n. 

The Glauber dynamics for the Ising model is a family of continuous-time 
Markov chains on the state space reversible with respect to the Gibbs 
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distribution, given by the generator 



(1.2) 



where is the configuration a with the spin at x flipped. The transition 
rates c(x, a) are chosen to satisfy finite range interactions, detailed balance, 
positivity and boundedness and translation invariance (see Section [2]) . Two 
notable examples for the transition rates are 

(i) Metropolis: c{x,a) = exp [2ha{x) +2(5<j{x)^y^^<j{y)^ A 1. 

(ii) Heat-hath: c(x, cr) = 1 + exp ^ — 2ha{x) — 2[5(j{x) Yly,^x ^iv) 

These chains have useful graphical interpretations: for instance, heat-bath 
Glauber dynamics is equivalent to updating the spins via i.i.d. rate-one 
Poisson clocks, each time resetting a spin according to the conditional dis- 
tribution given its neighbors. 

Ever since its introduction in 1925, the static properties of the Ising 
model and its Gibbs states, and more recently the Glauber dynamics for 
this model, have been the focus of intensive research. A series of break- 
through papers by Aizenman, Dobrushin, Holley, Shlosman, Stroock et al. 
(cf., e.g., [l][l5][l^|2l|[2§|25)[30}|32|[^^ starting from the late 1970's has 
developed the theory of the convergence rate of the Glauber dynamics to 
stationarity. It was shown by Aizenman and Holley [l] that the spectral 
gap of the dynamics on the infinite-volume lattice is uniformly bounded 
whenever the Dobrushin-Shlosman uniqueness condition holds. Stroock and 
Zegarlihski (35,37 38 proved that the logartihmic-Sobolev constant is uni- 
formly bounded provided given the Dobrushin-Shlosman mixing conditions 
(complete analyticity) . Finally, Martinelli and Olivieri [30,31 obtained this 
for cubes under the more general condition of strong spatial mixing. This 
in particular established O(logn) mixing throughout the uniqueness regime 
in two-dimensions. See the excellent surveys ^28,,29j for further details. 

To conclude this collection of seminal papers that altogether established 
0(log n) mixing throughout the regime of strong spatial mixing, it remains 
to pinpoint the asymptotics of the mixing time and determine whether or 
not there is cutoff in this regime. 

The cutoff phenomenon was first identified for random transpositions on 
the symmetric group in [12], and for the rifSe-shuffle and random walks on 
the hypercube in ^ . The term "cutoff" was coined by Aldous and Diaconis 
in [4], where cutoff was shown for the top-in-at-random card shuffling pro- 
cess. See [gIPz 34 and the references therein for more on the cutoff phenom- 



enon. In these examples, and most others where cutoff has been rigorously 
shown, the stationary distribution has many symmetries or is essentially 
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one-dimensional (e.g. uniform on the symmetric group \i2j, uniform on the 
hypercube ^ and one-dimensional birth-and-death chains |14|). Even for 
random walks on random regular graphs (where the stationary distribution 
is uniform), cutoff was only recently verified [27]. 

In the context of spin systems, cutoff was conjectured by Peres (see [22[ 
Conjecture 1]) for the Glauber dynamics on any sequence of transitive graphs 
on n vertices where its mixing time is O(logn). More specifically, it was 



conjectured in 22 that cutoff holds for the Ising model on any d-dimensional 



torus (Z/nZ)'^ in the high temperature regime; see also 23 Question 8], 
where the special case of d = 1 (Ising model on the cycle) was emphasized. 

However, so far the only spin-system where cutoff has been established at 
some inverse-temperature /? > is Glauber dynamics for the Ising model on 



the complete graph [13 22 . There, the magnetization (sum-of-spins) is in 
fact a one-dimensional Markov chain whose mixing and cutoff govern that 
of the entire dynamics. While this result motivates the conjecture on cutoff 
for the Ising model on lattices, its proof fails to provide insight for the latter 
setting since the complete graph has no geometry to consider. 

1.2. Cutoff for spin systems on the lattice. Our proof of cutoff for 
the Ising model on the d-dimensional lattice, as stated in Theorems T]|2 



in fact applies to a broad class of spin systems (essentially any monotone 
or anti-monotone system), underlying geometries and boundary conditions. 
We demonstrate this by establishing cutoff for the anti-ferromagnetic Ising 
model and the gas hard-core model (see, e.g., p3| for definitions). 

Theorem 5. The cutoff result given in Theorem for the ferromagnetic 
Ising model also holds for the anti-ferromagnetic Ising model. 

Theorem 6. Let d > 1. The following holds for the Glauber dynamics for 
the gas hard-core model on (Z/nZ)'^ with fugacity (3 and periodic boundary 
conditions. If (3 is such that there is strong spatial mixing, then the dynamics 
exhibits cutoff at ((i/2Aoo) log n with window of Oi\og\ogn) , where Aqo is the 
spectral gap of the dynamics on the infinite volume lattice. 

Our methods also establish cutoff on more general spin-systems such as 
the Potts model, provided that the temperature is sufficiently high. This is 
discussed further in the companion paper \2Q\. 

Furthermore, our results are also not confined to the underlying geome- 
try of the cubic lattice and in fact Theorem [2] holds for any non-amenable 
translation invariant lattice (e.g. triangular/hexagonal/ladder lattice etc.). 

Similarly, the choice of periodic boundary conditions is not a prerequisite 
for establishing cutoff, though it does enables us to determine the mixing 
time in terms of the infinite- volume spectral gap. For arbitrary boundary 
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Figure 2. Sparse geometry of the update support along 
a short time interval (mininum subset of sites determining 
the final configuration, conditioned on the update sequence), 
simulated on a 500 x 500 square lattice at /3 = 0.4. Color- 
map highlights the last time a spin belonged to the support 
(black being earliest and white being latest). 

conditions we can still establish the existence of cutoff and in various spe- 
cial cases (e.g. under the all-plus boundary conditions) we can pinpoint its 
location as stated by the following results of the companion paper [26] . 

Theorem 7 ([26]). Let d > 1 and consider the Glauber dynamics for the 
Ising model on {Ij/nLY with arbitrary boundary conditions. If there is strong 
spatial mixing then the dynamics exhibits cutoff. 

Theorem 8 ([26^). Consider Glauber dynamics for the Ising model on 
(Ij/n'L)'^ with all-plus boundary condition. If there is strong spatial mix- 
ing then the dynamics exhibits cutoff at (Aoo A 2Ae)^^logn, where Aoo,Ae 
are the spectral gaps of the dynamics on the infinite-volume lattice and on 
the halfplane with all-plus boundary condition respectively. 
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1.3. Main techniques. To prove total- variation cutoff on the cubic lattice 
of dimension we introduce a technique for bounding the L^-mixing of the 
dynamics on (Z/nZ)'* from below and above in terms of a quantity that 
measures the L^-mixing of the dynamics on smaller lattices (Z/rZ)'^ with 



r X log^ n (see the definition of rrit in (3.1 )). 



A key ingredient in this technique is an analysis of the geometry of the 
"update support": the minimum subset of sites whose initial spins, condi- 
tioned on the update sequence in a given time period, determine the final 
configuration. We show that even for relatively short update sequences, this 
minimum subset of sites is typically "sparse" , i.e. comprises remote clusters 
of small diameter. Namely, the clusters have diameter 0(log^ n) and the 
pairwise distances between distinct clusters are all at least of order log^ n 



(see Definition 3.3 of a sparse set). Fundamentally, this is ensured by the 
monotonicity of the system and the strong spatial mixing property, as the 
dynamics on cubes mixes in timescales much smaller than their diameter. 
This is illustrated in Figure [2] where the clusters are shown in white. 

We reduce the analysis of the L-'^-mixing of the dynamics on (Z/nZ)^ to its 
projection onto sparse supports. For such subsets, the speed of propagation 
of information between the distant clusters makes the projections onto them 
essentially independent. Using some additional arguments, from here we can 
translate the above to a product chain on smaller tori (Z/rZ)"^. Finally, the 
use of log-Sobolev inequalities provides a tight control over the L^-mixing 



in terms of the aforementioned L -quantity, as stated in Theorem 3.1 



2. Preliminaries 
Consider the Ising model on set of spins V, let n be its Gibbs distribution 



as given in ( 1.1 ) and a E O = {±1}^ be a configuration. For A C V, we will 
let (t(A) and fiA denote the projections of a and /U onto {±1}^ respectively. 
We let /i^Y denote the measure on A given the boundary condition r, that 
is, the conditional measure //a(' | ctoa = t)- 

2.1. Glauber dynamics for the Ising model. The Glauber dynamic for 
the Ising model on the lattice V = Z*^, whose generator is given in (1.2), 



accepts any choice of transition rates c{x, a) which satisfy the following: 

(1) Finite range interactions: For some fixed R > and any x E if 
o", o"' E agree on the ball of diameter R about x then c(x, a) = c{x, a'). 

(2) Detailed balance: For all a £ 0, and x £ V, 

^^=e^v{2ha{x)+2Mx)Y^a{y) 
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(3) Positivity and boundedness: The rates c{x,(t) are uniformly bounded 
from below and above by some fixed Ci, C2 > 0. 

(4) Translation invariance: If ex = cr'{- + i), where i £ V and addition is 
according to the lattice metric, then c(x, a) = c(x + i, a') for all x £V. 

The Glauber dyanmics generator with such rates defines a unique Markov 
process, reversible with respect to the Gibbs measure fiy. For simplicity of 
the exposition, the proof is given for the cases of heat-bath Glauber dynamics 
or Metropolis Glauber dynamics. The results for other choices of transition 
rates follow with minor adjustments to the arguments. 

We will let IPo-o(') s-^d IEcro[-] denote the probability and expectation con- 
ditioned on the Glauber dynamics having initial configuration ctq. 

2.2. Mixing, spectral gap and the logarithmic- Sobolev constant. 

The (total-variation) distance is perhaps the most fundamental notion of 
convergence in the theory of Markov chains. For two probability measures 
z^i, f2 on a finite space the total- variation distance is defined as 

\Wi - z^2||tv = max\ui{A) - U2{A)\ = - J^2{x)\ , 

Acq 2 ^ — ' 

i.e. half the L^-distance between the two measures. For an ergodic Markov 
chain (Yj) with stationary distribution z/, the mixing-time notion t^ix is 
defined with respect to ||P(lt £ ■) ~ /^IItv- 

The spectral gap and log-Sobolev constant of the continuous-time Glauber 



dynamics are given by the following Dirichlet form (see, e.g., 28 33]): 



^ = i^f.r^' «^=i^f/777Y' (2-1) 
/ Var(/) / Ent(/) 

where the infimum is over all nonconstant / G L^ifJ-) and 

<T,X 

Ent(/)=E [/2(a) log (/2(a)/E/2(cT))] . 



It is well known (see e.g. |3||10]) that for any finite ergodic reversible Markov 
chain < 2as < A and < tMix(l/e). In our case, since the sites are 
updated via rate-one independent Poisson clocks, we also have A < 1. 

By bounding the log-Sobolev constant one may obtain remarkably sharp 
upper bounds not only for the total-variation mixing-time but also for the 
L^-mixing (cf., e.g., [8-11,33). The following theorem of Diaconis and 



Saloff-Coste [To[ Theorem 3.7] (see also [33], [s] Chapter 8]) demonstrates 
this powerful method. 
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Theorem 2.1. Let (Yt) be a finite reversible continuous-time Markov chain 
with stationary distribution u. For any x with v{x) < and any s > 0, 

Wx{Xs G •) - ^L^iu) ^ exp (l- X(s- ^ log log 



UiX 



2.3. Strong spatial mixing and logarithmic-Sobolev inequalities. As 

noted in the introduction, bounds on the log-Sobolev constant of the Glauber 
dynamics for the Ising model were proved under a variety of increasingly gen- 
eral spatial mixing conditions. We will work under the assumption of strong 
spatial mixing (or regular complete analyticity) introduced by a Martinelli 
and Oliveri (SOl as it holds for the largest known range of (3. 



Definition 2.2. For a set A C Z'^ we say that SM(A,ci,C2) holds if there 
exist constants ci, C2 > such that for any A C A, 

< cie-^2distfe,A) ^ 



sup 



[ma)a - (^A )a 



where the supremum is over all y G dA and r G {±1}'^^ and where (/^]\^)a is 
the projection of the measure fi\ onto A. We say that strong spatial mixing 
holds for the Ising model with inverse temperature (3 and external field h on 
Ia'^ if there exist L, ci, C2 > such that SM(Q, ci, C2) holds for all cubes Q of 
side-length at least L. 

The above definition implies uniqueness of the Gibbs measure on the 
infinite lattice. Moreover, strong spatial mixing holds for all temperatures 
when d = 1 and for d = 2 it holds whenever /i 7^ or /5 < /?c. As discussed in 
the introduction, this condition further implies a uniform lower bound on the 
log-Sobolev constant of the Glauber dynamics on cubes under any boundary 



condition r (see [28 , 30 31] ) . We will make use of the next generalization of 
this result to periodic boundary conditions, i.e. the dynamics on the torus, 
obtained by following the original arguments as given in j28] with minor 
aherations (see also [5][l7)[25][29] ) . 

Theorem 2.3. Suppose that the inverse-temperature {3 and external field h 
are such that the Ising model on has strong spatial mixing. Then there 
exists a constant a* = a^{l3, h) > such that the Glauber dynamics for the 
Ising model on (Z/nZ)*^ with periodic boundary conditions has a log-Sobolev 
constant at least a* independent of n. 

3. Reducing 1} mixing to local mixing 

In this section, we show that the (total-variation) distance of the dy- 
namics on the lattice from stationarity is essentially determined by the 
distance from stationarity of a projection of this chain onto smaller boxes. 
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More precisely, consider the continuous-time Glauber dynamics (Xt) for 
the Ising model on Z^, the d-dimensional lattice with side- length n and 
periodic boundary conditions, and let denote its Gibbs distribution. Fur- 
ther, consider such a chain on a smaller lattice, namely (X^ ) on for 
r = 3 log^ n, and let /x* denote its stationary distribution. (We actually 
have a lot of freedom in the choice of r, e.g. any poly-logarithmic value 
which is at least log n for some fixed 6 > would do; this will prove 
useful later on for relating the cutoff location to the spectral gap of the 
dynamics on the infinite- volume lattice.) Within this smaller lattice we con- 
sider a d-dimensional box B with side-length 2 log^ n (the location of the 
box B within does not play a role as the boundary is periodic). Define 

xnt = max\\F.,,{XUB) G •) " /^sIL^f^M - (3-1) 

where X^^B) and /U^ are the projections of X^ and /x* resp. onto the box 
B. The following theorem demonstrates how the mixing measured by 
the quantity n'^xn^ governs the mixing of (Xt). 

Theorem 3.1. Let (Xt) be the continuous-time Glauber dynamics for the 



Ising model on Z^, and define xwt as in (3.1). The following then holds: 

1. Let s = s(n) and t = t{n) satisfy (lOd/a*) loglogn < s < log'^'^^n and 
< t < log^''^ n. For any sufficiently large n, 

uiax\\F,,{Xt+s G •) -/^IItv < \ (exp((n/log5n)'^m?) - l) + 3n-9'^ . 
In particular, if {n / log^ n)'^ as n ^ oo for the above s,t then 

limsupmax WF^^^Xt+s G •) " /^IItv = • 

2. // (n/ log^ n)°'mj oo for some t > (20(i/a*) loglogn, then 

liminf max ||P2;„(Xt G •) — A'lLv = 1 . 

Remark 3.2. It will be useful to apply Part 1 of the above theorem to 
lattices of varying sizes. Indeed, we will show that if (Xt) is the continuous- 
time Glauber dynamics for the Ising model on Z^ with 

log^ n < m < n , 



and with s,t,mt as in Theorem 3.1 (e.g., s,t < log^'^^n, the box B in the 
definition of mt given in (3.1) has side length 21og^n, etc.), then 



max \\F.,,{Xt+s S •) - /^IItv < I fexp (m'^m?) - l)'^' + Sn-^'^ . (3.2) 

Xo I \ ' J 
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The upper and lower bounds stated in the above theorem appear in Sub- 
sections [3]2] and |3]3] respectively. We begin by describing two key ingredients 
in the proof of the L^-Lp' reduction that enable us to eliminate long-range 
dependencies between spins and break down the lattice into smaller inde- 
pendent blocks. First we analyze which spins effectively influence the final 
configuration at some designated target time and characterize the geomet- 
ric structure of components comprising such spins. Second, we introduce 
the harrier- dynamics^ a variant of the Glauber dynamics that separates 
the lattice into weakly-dependent blocks by surrounding each one with a 
periodic-boundary barrier. 

3.1. Eliminating long-range dependencies. Consider some time point 
t = t{n) X logn just prior to mixing, and let s = s{n) x log log n be a short 
time-frame. Our goal in this section is to bounds the L^-distance of the 
Glauber dynamics from equilibrium at time t-|-s in terms of the L^-distance 
projected onto sparse subsets of the spins. 

Definition 3.3 {Sparse set). Let log^n < m < n. We say that the set 
A C is sparse if for some L < (n/log^n)*^ it can be partitioned into 
components Ai, . . . , such that 

1. Every Ai has diameter at most log^ n in Z^. 

2. The distance in Zj^ between any distinct Ai,Aj is at least 2dlog^ n. 
Let S = S{m) = {A C ZJ^ : A is sparse}. 

Theorem 3.4. For log^ n < m < n let (Xt) be the Glauber dynamics on 
Z^ and n its stationary measure. Let (lOd/a*) log logn < s < log^^^n and 
t > 0. Then there exists some distribution u on S{m) such that 

||P,o(Xt+, G •) - /"IItv < ^ Wx,XXt{^) G •) - /"aIItv M^) + 3n-'°^ . 

Proof. To prove the above theorem, we introduce the following variant of 
the Glauber dynamics which breaks down the dynamics into smaller blocks 
over which we have better control. 

Definition 3.5 {Barrier- dynamics). Let (Xt) be the Glauber dynamics for 
the Ising model on Z^. Define the corresponding barrier- dynamics as the 
following coupled Markov chain: 

(1) Partition the lattice into disjoint d-dimensional boxes (or blocks), 
where each side- length is either log^ n or log^ n — 1. 

(2) For each such box B, let B~^ be the d-dimensional box centered at B 
with side-lengths log^ n -\- 21og^''^ n, e.g., if B has side-length log^ n 

5+ = y {u : \\u - v||oo < log^/^ n} . 
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Let V's be a graph isomorphism mapping onto some block C"*" 
(and B onto C C C"*"), where the C"^ blocks are pairwise disjoint. 
(3) Impose a periodic boundary condition on each (to be thought 
of as a barrier surrounding it), and run the following dynamics: As 
usual, each site u in ZJ^ receives updates according to a unit rate 
Poisson clock. Updating u at time t via a variable / ~ U[0, 1] in the 
standard dynamics implies updating every v = tpBiu) (for some B 
with u G B+) via the same update variable I. 

The above Markov chain gives rise to the following randomized operator 
Qs (for s > 0) on {±1}^™. Given an initial configuration xq for Z^, we 
translate it to a configuration for the blocks in the obvious manner, and 
then run the barrier-dynamics for time s. The output of the operator Qg 
is obtained by assigning each u G Z^ the value of i/jb{u), where B is the 
(unique) block that contains u. In other words, we pull-back the configura- 
tion from the Cs onto ZJ^ (while discarding the spins of the overlaps). 

To simplify the notations, we identify the blocks C, with B, B~^ when- 
ever there is no danger of confusion. 

Note that for the mixing-time analysis, we are only interested in the 
behavior of the dynamics up to time O(logn), and the above parameters 
were chosen accordingly. Indeed, the next lemma shows that the barrier- 
dynamics can be coupled to the original one up to time (logn)^/^ except 
with a negligible error-probability. 

Lemma 3.6. Let to = (logn)"^/^. The barrier- dynamics and the original 
Glauber dynamics are coupled up to time to except with probability n~^^'^. 
That is, except with probability n~^^'^, for any Xo we have Xg = Gs{^o) 
simultaneously for all s < to- 

Proof. Let {Xtj denote the Glauber dynamics on Z^ and let {Xt) denote the 
barrier-dynamics with corresponding blocks Bi and B^ . Apply the afore- 
mentioned coupling between the two processes, where the original Glauber 
dynamics runs as usual, and the barrier-dynamics uses the same updates for 
each of its sites. That is, if site u is updated in the original dynamics via 
a uniform real / ~ C7[0, 1], we update it in every that contains it using 
the same / at the same time. 

Consider some box B and its block B^. Clearly, the barrier-dynamics on 
is identical to the original dynamics until it needs to update the 
boundary of B~^ (in which case (Xf) has periodic conditions whereas (Xf) 
uses external sites in the lattice). 

Therefore, a necessary condition to have Xt{v) / Xt{v) for some v & B 
is the existence of a path of adjacent sites ui,. . . ,ue connecting v to dB'^, 
and a sequence of times ti < . . . < t^ < t such that site Uj was updated 
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at time ti (note that £ > log'^/^n by definition). Summing over all {2dY 
possible paths originating from v, and accounting for the probability that 
the i corresponding rate 1 Poisson clocks fire sequentially before time t < to 
(while recalling that Iq = o{i)) it then follows that 

J^P (Ui {Xt{B,) / Xt{Bi)^^ < ton'^ Yl (2d)^P(Po(to) > i) 

t<to e>log^^^ n 

where the inequalities hold for any large n (with room to spare). ■ 

In light of the above lemma, we can focus on the barrier-dynamics for 
the sake of proving Theorem |3.4[ Crucially, suitably distant sites evolve 
independently in this new setting. 

The random operator Qg is determined by the random update sequence Wg 
(each update is a tuple {uj,tj,Ij), where uj is the site that was updated, tj is 
the time of update and Ij is the unit variable determining the update result). 
In other words, for any such sequence Ws there exists some deterministic 
function gw^ : {±1}^™ — > {±1}^'" so that ^s(x) = gw_^{x) for all x. Further 
note that is monotone, by the monotonicity of the Ising model. We use 
the abbreviated form P(VFs) for the probability of encountering the specific 
update sequence Ws between times (0,s). 

Definition 3.7 {Update support). Let Wg be an update sequence for the 
barrier-dynamics between times (0,s). The support of Wg is the minimum 
subset C such that Qs{x) is a function of x{Aws) for any x, i.e., 

9Ws{x) = fws{x{^wj) for some fw, : {±1}^^^ ^ {±1}< and all x. 

In other words, v ^ A.]^^ if and only if for every initial configuration x, 
modifying the spin at v does not affect the configuration gwsi^)- This 
definition uniquely defines the support of Wg. 

Using this notion of the support of updates in the barrier-dynamics, we 
can now infer the following upper bound on the L^-distance of the original 
dynamics to stationarity. 

Lemma 3.8. Let (Xt) be the Glauber dynamics on ZJ^, and Wg be the 
random update sequence for the barrier-dynamics along an interval (0, s) 
for some s < log^*^^ n. For any xq and t > 0, 



lOd 
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Proof. Let Xt be the Glauber dynamics at time t started from Xq = xq, 
as usual let = {±1}^™ and let Y G be distributed according to /i. 
Recalling that gw^ denotes the deterministic function associated with an 
update sequence Wg for the barrier-dynamics in the interval (0,s), for any 
random configuration X € fim we have 

\\F{g,{X) G •) - HGsiY) G OIItv = max mGs{X) G A) - P(g,(y) G A)] 



max / mgwAX) G A) -F{gwAY) G A)] dF{Ws 



ACftr, 

Since gw^i^) = fWsi^i^Ws)) by definition of A^^, the above is at most 



max [F{fwM{AwJ) G A) - P(/^^(y(AH/J) G A)] dF{Ws) 

AG' £m 

< J \\F{X{AwJ G •) -P(nAwJ G OILv 

= 1 \\F{X{AwJ G ■) - /^^A,,, II,, c?P(W^s) , 

where the inequality in the second line used the fact that when taking a 
projection of two measures thei r tot al-variation can only decrease. 

Since s < log^''^ n, by Lemma 3.6 we can couple Xt^g and Qs{Xt) together 
except with probability n~^^'^ and hence 



\\F{Xt+sG-)-Figs{Xt)€-)\\^^<n 



-lOd 



Similarly, by Lemma 3.6 we can couple Gs{y) with the Glauber dynamics 
run from Y for time s (having the stationary distribution fi) and hence 



lOd 

TV ^ " 



||p(g,(y)G-)-/^llTv<" 

Combining these estimates, it follows that 
\\F{Xt+g e.)-f,\\^^.< mOgiXt) G •) - FiGsiY) G •) IItv + 2n''°'^ 

< I \\F{XtiAwJ G •) - MA^, dFiWs) + 2n-io^ , 
as required. ■ 

Thus far, we have established an upper bound on the L^-distance between 
Xt+s and in terms of Aiy^ C Z^, the support of the update sequence in 
the barrier-dynamics operator Gs along the interval (0, s). We now wish to 
investigate the geometry of the set of sites comprising Aws for a typical 
update sequence Ws- The following lemma estimates the probability that 
Aw^ is sparse, as characterized in Definition 3.3, Figure [3] shows a realization 
of the update support becoming sparser with time. 
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Figure 3. Evolution of the update support over increasing 
time intervals, simulated on a 500 x 500 square lattice at 
inverse-temperature [5 = 0.4. Highlighted regions correspond 
to the components comprising the sparse support. 



Lemma 3.9. Let Qg he the harrier- dynamics operator on TU^^, let Wg he the 
update sequence up to time s for some s > (lOd/os) loglogn, and S he the 
collection of sparse sets of Z^. Then F{Aws ^ S) > 1 — n~^^'^ for any 
sufficiently large n. 

Proof. First, consider a single block B in the barrier-dynamics on Zj^ and 
for simplicity let vr stand for ^ib+i the stationary distribution on . Let 
Eb denote the event that Aiy^, D B ^ (}> for a random update sequence Wg. 

Observe that, by definition, the following holds for all s: li B ^ B are 
two distinct blocks, then modifying the value of a spin u ^ B m. the initial 
configuration Xq can only affect Xs{B^) provided that B^B are adjacent 
(otherwise the barrier around B prevents the effect of this change). Hence, 
when assessing whether u € it suffices to consider the projection of Xg 
onto the block B and all of its neighboring blocks; let N{B) denote this set 
of S'^ blocks. 

Let {X^) and {X^) be two instances of the barrier-dynamics restricted to 
B'^ starting from the all-plus and all-minus states respectively, coupled via 
the monotone coupling (note that the restriction to one block turns these 
into the standard Glauber dynamics). Clearly, if these two chains coalesce 
under the update sequence Wg at some point < t < s then either one of 
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them predicts Xs{B^) regardless of the value of Xq{B^), and so 

^BeN{B) ^ 
where the last inequality is by symmetry. As the system is monotone, 

ueB+ 

(3.3) 

Moreover, for any u £ 

[xtiu] / x-(n)) < ||P(X+ G •) - vtIItv + ||P(^7 e •) - vtIItv 

< i||P(X+ G •) - vtIIl^w + ^||P(X- G •) - vtIIl^w , 



and by Theorem 2.1 if the all- plus state 1 has stationary measure at most 
e~^ (clearly the case for large n) then for any s > 



e ■) - IT 



< exp ( 1 — A ( s — - — log log 



where A and as are the spectral gap and log-Sobolev constant resp. of the 
Glauber dynamics on B+. Recalling that 7r(l) > = 2-(i+"(i))i°g"'" 

and that X > as > a*, the assumption on s gives that 

1 8d 
s > — loglog(l/7r(l)) + y log log n 



for any sufficiently large n, and in this case 

2W 



||P(X+G-)-vrk2M<3(logn)-8'=', 



By the exact same argument we have 

||P(X7 gO-vtIIl^w <3(logn)-8^ 
and it now follows that 

P(x+(n)/X7(n))) <(3 + o(l))(logn)-6'^<4(logn)-6^ 

u£B+ 



where the last inequality holds for large n. Combining this with (3.3) yields 
that for any large n, the following holds with room to spare: 

P(i^B) <(ilogn)-6^. 

This estimate will now readily imply a bound on the number of com- 
ponents in the support of Wg. Let denote the following event: There 
exists a collection B oi L> (n/log^n)'^ blocks, such that Eb holds (that is, 
Ci B ^ 9) for all B € B, and the pairwise distances in blocks between 
the blocks in B are all at least 4. We claim that P(£''') < n^^^'^ for large n. 
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To see this, first notice that if a block B has a distance of at least 4 
blocks from a set of blocks B (i.e., for any B' £ B', no two blocks in 
N(B) and N(B') are adjacent), then the variable Ie^ is independent of 
{1^^, : B £ B}. Indeed, the initial configuration on B can only affect the 
outcome of Xs{N(B)), and by definition this outcome is derived from the 
initial configuration via the updates of the sites Ub = U^g^^^^B"*". Our 
assumption on the distance between B and B precisely implies that the 
above Ub is disjoint to any Ub' for B' ^ B in B, and the statement now 
follows from the independence of updates to distinct sites. Hence, as the 
total number of blocks is (1 + o(l))(n/log^ n)'^, for large n we have 

Now suppose that there is a sequence of blocks, (Sj,, , Bi^ , Bi^ , ■ ■ ■ , Bi^) for 
some ^ > ^0 = log such that for all k, 

(1) The distance in blocks between Bii__-^,Bii, is at most 3d (i.e., they 
are the endpoints of a path of at most 3d+ 1 adjacent blocks). 

(2) We have Aw, n / 0. 

Clearly, if the distance in blocks between some B,B' is at least 3d+l, then 
since every block has side-length at least log^ n — 1, the distance between 
any two sites u £ B and v £ B' is at least 2d\o^ n for sufficiently large n. 

Observe that if Ap^/^ ^ 5, then either the event holds or a sequence 
of blocks as described above must exist. Indeed, consider some A^y^ that is 
not sparse, and partition it into components as follows: 

u G S n Avi/^ and u' B' f] Ayy, belong to the same component 
<^=^ the distance in blocks between B, B' is at most 3d. 

The number of components is clearly at most (n/log^n)'^, otherwise the 
event holds. Furthermore, by definition (as argued above), the distance 
between any two distinct components is at least Idlog^ n. Hence, the as- 
sumption that Ai4/^ ^ S implies that some component Ai must have a diam- 
eter larger than log^ n. In particular, there are two sites u, u' G Ai belonging 
to B,B' respectively, such that B,B' have distance of at least logn blocks 
between them. Moreover, by the way we defined the component Ai there are 
blocks B = Bj^, . . . , Bj^ = B' such that Bj. n ^ and the distance in 
blocks between Bj. , Bj-_^-^ is at most 3d. As the assumption on u, v! implies 
that I > j^logn, this sequence satisfies the required properties. 

It therefore remains to show that, except with probability n~'^^'^, there 
does not exist a sequence of blocks satisfying the above properties ([T|,([2]). 

The above described sequence in particular contains a set B of at least 
£/7'^ blocks, whose pairwise distances are all at least 4, and A\y^ r\ B ^ % 
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for all B £ B (for instance, take B = {-Bj^}, process its blocks sequentially, 
and for each B delete from B any B' ^ B whose distance from B is less 
than 4, a total of at most 7^^ — 1 blocks). By the above discussion on the 
independence of the events Eb for such blocks, as well as our estimate on 
¥(Eb), we deduce that for any large n 

F{Aw, n 5 / for all 5 G ^) < log n) "^"^^^^^ . 

Clearly, there are at most n^{6d+l)'^^ sequences of blocks {Bi^, Bi-^, . . . , 
where any two consecutive blocks are of distance (in blocks) at most 3d. 
Hence, the probability that there exists a sequence with the aforementioned 
properties ([l]),([2]) is at most 

(6d+ 1)^^(1 log n)-*^j < n'^(logn)~(5'^/^')^o < ^'^-7"^ log log n 

i>io 

for any large n, as required. ■ 



Theorem 3.4 is now obtained as an immediate corollary of Lemma 3.8 and 



Lemma 3.9, with u given by u{A) = F{{Ws : Aw, = A}) for AeS. 



3.2. Proof of Theorem 3.l[ Part 1: upper bound on the distance. 



Let A C be a sparse set and U^^^Ai be its partition to components as 
per Definition 3.3, where L < m'^ A (n/log^n)'^. Letting dist(-,-) denote 



the distance according to the lattice metric of 'Z'^, put 

Af = {v : dist(^i, v) < log^/^ n} , 

and (recalling that the diameter of Ai is at most log^ n) let be an isometry 
of Af into a box Bi of side-length 2 log'^ n inside Z^, the torus of side-length 
r = 31og^n. A crucial point to notice is that the sets {^^} are pairwise 
disjoint, since the distance between distinct Ai,Aj is at least 2dlog^n. 

Let {Xf ) denote the product chain of the Glauber dynamics run on each of 
these L copies of ZJ? independently, and let /i* denote its stationary measure. 
Given some initial configuration xq on Z^, we construct Xq for {X^) in 
the obvious way: For each v G Af , assign XQ{Tpi{v)) = xo{v), while the 
configuration of each Zj? \ tlJi{Af) can be chosen arbitrarily (e.g., all-plus). 

Let t < log'^'^^n. We claim that we can now couple (Xg) with {X*) such 
that Xs{A) = X* (UipiiAi)) for all < s < t except with probability n'^^'^. 



Indeed, by repeating the argument of Lemma 3.6 if we let {X*) use the 
same updates of {Xs) on UAf and run independent updates elsewhere, then 
each Ai can be coupled to ipi{Ai) on [0, t] except with probability n~^^'^ (here 
too the distance between Ai and dAf is at least log^''^ n). Since the subsets 
{Af} are disjoint, distinct V'i(^i) indeed obtain independent updates in this 
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manner. Summing the above error probability over the L copies we deduce 
that the coupling exists except with probability Ln^^^'^ < n^^^'^. 

We claim that by the strong spatial mixing property, for any large n 

ll/^A - ^(uv,(A,))IItv < n-^"^ (3.4) 

after one translates the sites in the obvious manner (that is, when comparing 
the two measures we apply the isometries ipi on the inputs). 

It is easy to infer the above statement under the assumption that there is 
strong spatial mixing SM(A, ci, C2) for any subset A (so-called complete ana- 
lyticity, as opposed to our assumption that strong spatial mixing SM((5, ci, C2) 
holds for sufficiently large cubes). Indeed, recall that for some fixed ci, C2 > 
we have — /^X^IItv < ci exp(— C2 dist(A, y)) for all r G {±1}^^ and 
y £ dA. By taking the boundary conditions ri and T2 to be distributed 
according to /i and /u* about UdAf and U^^pi{Af) respectively, we can 
transform ri into T2 via a series of at most m'^ spin flips. Each flip amounts 
to an error of at most ci exp(— C2 log^''^ n), leading to (3.4). 



Since we only have regular complete analyticity at hand, we will obtain 
(3.4) via monotonicity and log-Sobolev inequalities. Take t = (logn)^/'^ and 



apply Theorems |2.1| and 2.3 to get 



{Xt G •) - /^IItv < W^xoiXt G •) - ^^\\LH^,) 



< exp ^1 - A - ^ loglog;U*(2;o) 

where A, Us are the spectral gap and log-Sobolev constant of (Xt) respec- 
tively. Since loglog(l//i*(o")) > {d + o(l))logn for all a and since it holds 
that A > Os > CKg > 0, for a suitably large n we have 

max \\¥^^{Xt G •) " /^IItv < ^"^""^ • (3-5) 

XQ 



It is well known (see, e.g., |33 Lemma 2.2.11] and also [s", Chapter 8]) 



that the spectral gap (respectively log-Sobolev constant) of a product chain 
is equal to the minimum of the spectral gaps (log-Sobolev constants). In 
particular, the log-Sobolev constant of (X^) is at least a*, hence similarly 



max||P,.*(X;E.)-/"l,v^^"'''- (3-6) 
Therefore, if y is a configuration on A distributed according to /xa then 



+ ||P,„(V(Xi(A)) G ■)-F^^^{X:{uMA^)) G •)! 



+ 



^,2(X;(UV;i(A)) G-)-^(U^,(A,)) 



< n 

TV 



TV 

9d 
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where the last inequaUty holds for large n, combining (3.10), (3.5) and (3.6) 
with the fact that projections can only reduce the total variation distance. 
Altogether, letting T = UBi and abbreviating Fx*{X^ G •) by tt*, we have 



(3.7) 



We next seek an upper bound for the last expression, uniformly over the 
initial configuration Xq. 



bu-,bL 



< - E 



1 



2 \ 1/2 ^ 

~ 2 



VTr 



1 



2 

1/2 



VTr 



1 



(3., 



where the inequality was by Cauchy-Schwartz, and the last equality is due to 
the fact that E^* [TTr/z^r] ~ ^- Now, since (X*) is a product of L independent 
instances of Glauber dynamics on ZJ?, we infer that 



E 



Mr 



i=l ^B. i=l 



+ 1 , 



and recalling the definition of rrif in (3.1), it follows that E^* Kr//^rl 
most (m| + 1)^. Plugging this in (3.8), 



m^x ||p,.„ (x;(r) G •) - IItv < 2 (K + 1)^-1 

<:^(exp(Lm?)-l) 



1/2 



1/2 



Altogether, using (3.7), we conclude that for any A G 5 and t < log^''^ n, 
m^ax||P,„ {Xt{A) G •)- /^a|Itv < ^ (exp (Lm^) - l)'/' + 271-9-^ . 



At this point. Theorem 3.4 implies that for any s,t with < t < log^^^n 
and (lOd/a*) log log n < s < log^^^ n and for every large n. 



x'o(^t+s G •) - /^IItv < ^ (exp (Lmf) - l)^^^ + 3n 



-9d 



Recalling that L < m'^ A (n/log^n)'^ concludes Part 1 of Theorem 



3.1 



(3.9) 

and 



also establishes the inequality in Remark 3.2 
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3.3. Proof of Theorem 3.1 , Part 2: lower bound on the distance. 

Let 

r = 31og^n, L=[n/r\'^, 

and let Ai , . . . , C be a collection of d-dimensional boxes of side-length 
|r satisfying ||u — f ||oo > r/?, for any u ^ Ai and v S Aj with i ^ j. Similar 
to the notation of the previous subsection, for each i G {1, . . . , L} we define 



Af = {v : dist(^i,t;) < r/6}. 

Denote the unions of these boxes by A = uf^^^Ai and 

Let , . . . , be a sequence of disconnected d-dimensional boxes of 
side-length r and let F"*" denote the graph of their union. Let V'i be an 
isometry mapping Af to Bf and let ^p be the isometry that maps A"*" to r"*" 
such that its restriction to any individual Af is ipi. We define Bi = il^{Ai) 
and r = U^]^i?i. For a configuration X on some A' C A'^ we let will denote 
i^{X) as the corresponding configuration on ip{A'). 

We couple the Glauber dynamics on ZJ^ and as follows: Whenever 
a site u G A"*" receives an update via some unit variable / we also update 
the site ip{u) using the same / and a periodic boundary condition on its 
corresponding box Bf. Denote the dynamics induced on as X^, and let 
/i* be its stationary distribution. The above defined coupling satisfies the 
following claim. 

Claim 3.10. Let {Xt) and {X*) be the above coupled Glauber dynamics 
on and F+ respectively. Suppose Xq,Xq satisfy ilj{XQ{i^'^)) = Xq(F+). 
Then with probability at least 1 — n~^^'^, for all < s < (logn)^^^ we have 

iPiXsiA)) = X:(F) . (3.10) 



Furthermore, 



max||P^.„(X, G 

XO 



— /^IItv — I 



,{X:{T)e-)-fi*r\\^^-4n 



Proof. Equation ( 3.10 ) holds by a simple adaption of the proof of Lemma 3.6 



as the initial conditions and updates agree on Af and Bf for all i, and since 
each Bi is distance r /6 = ^ log^ n from the boundary of Bf . 

For the second statement of the claim, we repeat the argument that 
yielded inequality (3.4) in the previous subsection. Take t = (logn)^/^ and 



apply Theorems 2.1 and 2.3 to get 



■ 2^0 



{Xt G 



HI 



< exp 1 - A t 



4a, 



log log fi*{xo] 
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where A and as are the spectral gap and log-Sobolev constant of (Xt) re- 
spectively. This yields the following for any sufficiently large n: 



max I 

XO 



XO 



< n 



(3.11) 



Since the log-Sobolev constant of the dynamics on is at least a*, 

max||P,*(X; G < n'^^". (3.12) 

Now, if y G {±1}^ is distributed according to /xa then 

||P(V^(y) G •) -/^fIItv < \m{y) e ■) -Fxoii^ixm) G OIItv 
+ ||p,,(v^(Xt(A)) G •) -iP..5(^t*(r) G Oil,, 

+ ||P,*(X^^)G•)-;"^||,v^3n-lO^ 

where the last inequality added (3.10), ( 3.11| ) and (3.12) to the fact that 
projections can only reduce the total variation distance. It follows that 



max IIP,, (X, G •) 



^^\\^^ > max||P,„(V^(X,(A)) G •) " ^ O 



XQ 



>max||P,*(X:(r)G-)-/^^| 



where we used the facts that ||P(V'(y) G •) — ^J^y-Wtv 
||P(V^(X,(A)) G •)-¥,*{Xl{T) G OIItv < n-iO"^. 



< 3n 



An 



-Wd 



and that 



Recall that the Bi boxes have side-length 



2 log n, matching the 



boxes B in the definition (3.1) of m^. Let Xg = XQ{t) be a configuration on 
a box B of side-length 31og'^n which achieves rrif, i.e. 

mt=||P,2(x;(i?)G-)-/UB|' 
We define i.i.d. random variables 



y, ^ ^[Xl{B,) = V,\Xl{Bt)=xt) 

where the C/j are i.i.d. configurations on Bi distributed according to 
As the dynamics on different tori are independent it follows that the Yi are 
independent. As we will soon show, these random variables provide crucial 
insight into the mixing of [Xt) in the L^-distance. First we need to obtain 
some estimates on their moments. Clearly, 

[Xl{B^i) = h\Xl{Bt) = 



and 



Vary 



E 



{X*{Bi) = bi\X*{Bt) = xl) 



tJ-B^bi) 



m 



t ■ 



24 



EYAL LUBETZKY AND ALLAN SLY 



Moreover, by a standard L°° to reduction (cf., e.g., 33 ), 

\\Yi - llloo = IIP G • I X*(Bt) = xl) - /i^,J|Loo(^^ ) 



<||p(x;/2(i3^)G-|^o*(^^ 



2.1 



and 



and hence by Theorems 
absolute constant c > 0, yielding 



2.3 this is at most e ^ log log n £qj. gg^g 



E|1^-1P < ||K;-l|UVarl^ <e- 



!logn 2 



o(m?) . 



Define = log Yi . Taking Taylor series expansions gives 
EZi = E(Yi - 1) - ^E(y, - 1)2 + O (El^i - 1|3) = ■ 
and similarly. 



l-o(l) ^2 
2^t > 



E{Y, - 1)2 + O (ElFi - if) = (1 + o(l))m: 



it is 



We now derive the required lower bound on the L^-distance of (Xf) from 
stationarity at time t > (20d/a*) log log n, provided that lim„_^oo Lm.^ = oo. 
Consider a starting configuration Xq such that ip{Xo{Af)) = Xq for all i, 
and similarly choose Xq so that Xq{B^) = Xq for all i. By Claim 3 
sufficient to show that under our hypothesis 

||p,„(x;(r)G.)-/^rllTv^i- 

By the definition 1^'s, 

||p(x;(r)G.)-MrllTv 

= 1 \^.oiX:iBi) = h,...,x:{BL) = bL)-^*r{bi,...,bL)\ 



E n 

fei,...,6i 1=1 

L 



XQ 



{Xt{B,) 



i=l 



-E 



1=1 



l{Yi-l =E exp^Z, -1 



i=l 



where |a|~ denotes max{— a,0}. The random variables Zi are independent, 
||-^i||oo = o(l) and X^i^iVarZj — > oo, hence the Central Limit Theorem 
implies that 

lim^ti<^^=Af(0,l). 
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Since LEZi/y^Lm^ = — ^^^^ \J Lmf — oo, it then follows that YL\=\'^i 
converges in probability to — oo and therefore 



E 



exp 



i=l 



1, 



as required. 



4. Cutoff for the Ising model 
In this section we prove the main results, Theorems [T] and |2] We first 



describe how the L^-Lp' reduction from the previous section (Theorem 3.1) 
establishes the existence of cutoff. A refined analysis of the distance to 
stationarity (captured by the quantity mt in the aforementioned theorem) 
via log-Sobolev inequalities then allows us to pinpoint the precise location 
of the cutoff in terms of spectral gaps of the dynamics on tori of prescribed 
sizes. Finally, by applying this tool on tori of varying sizes we obtain as a 
biproduct that, as the tori side-length tends to infinity, these spectral gaps 
tend to Aoo, the spectral gap of the dynamics on the infinite- volume lattice. 
In turn, this yields the asymptotics of the mixing time in terms of Aoo and 
establishes Theorem IH 



4.1. Existence of cutoff. Theorem 3.1 already establishes cutoff for the 



Ising model on (Z/nZ)"^, although not its precise location. To see this, recall 

lose 

log=^ 



the definition of mt given in (3.1) and choose t* as follows 



t* =mi{t: mi < 



As before, the log-Sobolev inequalities of Theorems 2.1 and 2.3 imply that 
t* = O(logn). Let s = (lOd/a*) log log n. Since mt is a continuous funct ion, 
we have (n/ log^ n)'^m^, = log^^^'^ n = o(l), and so by Part 1 of Theorem 3.1 



m^ax||P^„(Xf.+, G •) - /^IItv < \ (exp ( log^-^d _ ^^^^ g^-M ^ ^ 

Next, the results of |18 imply that the mixing time of the Glauber 
dynamics for the Ising model on {'L/n'L)'^ has order at least logn, hence 
(by the above inequality) t*{n) is also of order at least logn. In particular, 
t* > (20d/as) log logn for any sufficiently large n, and since by definition 
(n/ log'^ n)'^mt. = logn, it follows from Part 2 of Theorem 



3.1 



that 



max ||P^.o(Xt* G •) - /^IItv = 1 " "(l) 

Xo 



This establishes cutoff at t* with a window of O(loglogn). 
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4.2. Cutoff location (asymptotics of the mixing time). To obtain the 
asymptotics of the mixing time, it remains to estimate the parameter t* 
introduced above, that is, to understand the threshold t{n) for n'^mf to tend 
to infinity faster than some poly-logarithmic function of n. 

In what follows, let A(r),as(r) be the spectral gap and log-Sobolev con- 
stant of the Glauber dynamics on a d-dimensional torus of side-length r. 



Lemma 4.1. Set cq = ^ , let log log n <t< log^/^n and r = 31og^ 
For n sufficiently large, 



e 



-A{r)t-cologlogn _ ^-9d ^ ^ g- A(r)t+co log log n _ ^^ -^-^ 



Proof. Let denote the Glauber dynamics on with periodic boundary 
conditions, let /x* be its stationary distribution and let Q* denote its state 
space. Since loglog(l/;U*((T)) > (Sd-I- o(l)) loglogn for all o" G $7* and since 



A(r) < 1 (vertices are updated at rate 1), another application of Theorem 2.1 
implies that for large n 

mt < maxexp ( 1 - A(r) [t - ^ loglog (l//i;(o-)) 

< g-A(r)t+2^?gploglogn. ^ g- A(r)t+co log log n ^ ^4^2) 

This establishes the upper bound on mt- Further, as t > {20d/a*X*) log log n 
it follows that r'^/'^mt < log'^'^n = o(l). 

A standard lower bound on the total variation distance in terms of the 
spectral gap (cf. its discrete-time analogue (23| equation (12.13)]) gives that 



g-A(r)t < 2 ||p(x; G •) - /^IItv for alH > . 
Set s = (lOd/Og) log log n. Applying Part 1 of Theorem 



3.1 to Z??f with these 



s,t (recalling Remark |3.2| and plugging in m = r in (3.2 )) gives 
g-AM(t+.)<2||p(x,%e.)-/"*Lv 

1/2 

< (exp (^r'^mf^ - l) + 6n~^'^ < 2r'^/'^mt + Gn'^'^ , (4.3) 



where the last inequality used the fact that for x < 1 we have e^ — 1 < 2a; 



and r'^mf = o(l). Rearranging equation (4.3) we have that 



mt > Q-MrKt+s)-log{2r'i/^) _ "^^^^^ > g- A(r)t-co log log n _ ^-9d ^^ ^^ 



Combining equations (4.2) and (4.3) completes the proof. ■ 

The following theorem now establishes the position of the mixing time in 
terms of A(r) with a window of O(loglogn). 
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Theorem 4.2. Let (Xt) the Glauber dynamics on the (Z/nZ)"', and set 

d 



t* 



t* 



15d 
atX* 



log log n, 



The following then holds: 



max 

Xo 



max 

XQ 



t +— -loglogn. 



1, 

0. 



Proof. We begin by applying Lemma 4.1 to t„ . The left-hand-side of (4.1) 
establishes that 



(n/log^ n)'^m^_ > (n/log'^n)'* ( n '*(logn)" 



-9d 



> (l + o(l))log3'^n, 



and in particular (n/log^ n)'^m^_ ^ co, hence by Part 2 of Theorem 



3.1 



max 



1 



Using the right-hand-side of ( |4.1[ ) for t = t+ — (lOd/a*) loglogn gives 
{n/log^nymj < (n/ log^ n)"^ • n-'^(log n)"^ < log""^ n = o(l) . 



Therefore, by Part 1 of Theorem 3.1 at time t^ = t + s we have 



max 



xo 



as required. 



4.3. Limit of spectral gaps. Theorem |4.2| established the location of the 
mixing time in terms of A(r), the spectral gap of the Glauber dynamics on 
the d-dimensional torus of side-length r. As commented at the beginning of 
Section |3] we are allowed some latitude in our choice of r. This provides a 
way of relating the spectral gaps for different values of r, ultimately proving 
that they converge to Aqo, the spectral gap on the infinite- volume lattice. 

Lemma 4.3. Let /? > 0, and let A(r) be the spectral gap of continuous- 
time Glauber dynamics for the Ising model on at inverse-temperature [5. 
If there is strong spatial mixing at inverse temperature (3 then there exists 
some A > such that 



|A(r) - A| < r 



-l/2+o(l) 



(4.5) 
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Proof. As noted in the beginning of Section [3j there is a lot of freedom in 
the choice of r = Slogan for the definition of rrij in (3.1), and the proofs 
hold as is (while resulting in slightly different absolute constants, e.g., cq in 
Theorem 4.2 etc.) for r = log^'^'^n with an arbitrary fixed 6 > 0. With this 



^etc.) for r 

in mind, fix some small 6 > and take ri = log^"*"^ n and ri < r2 < . 
In Theorem |4.2| we provided upper and lower bounds on the quantity 
IJ-W^y in terms of X{r),a*,X* and some absolute con- 

C{5) > such that for 



maXj.o ||Pxo(-'^t G • 

stants. As argued above, there exists some C 
every ri < r < rf the statement of the theorem holds with parameters 

d 



t* 



Cd 
atX* 



t*{n) 



log log n . 



2A(r 



log n . 



t* + 



Cd 
atX* 



log log n . 



Applying this theorem both on ri and r2 we must have t^{ri) > (r2) for 
sufficiently large n, and so 



d 



2A(ri) 



logn + 



Cd 
atX* 



log log n > 



d 



2X{r,: 



log n 



Cd 
a*X* 



log log n . 



Rearranging we obtain that 
A(ri 



Xir,) < 4g A(n)A(r2)loglogn ^ ^ y,^s 
a* A* logn 



where the last inequality holds for any sufficiently large n. As we can clearly 
reverse the role of ri and r2 it follows that for any large n. 



max |A(ri) — A(r2)| < 1/^+^ ^ 

ri<r2<rl 

By the above inequality, if r is a sufficiently large integer then 



Eh 

i=0 



A(r 



2i+lN 



< > r 

1=0 



oo . 



(4.6) 



(4.7) 



Combining equations (4.6) and (4.7) establishes that {A(r)}^o converges to 
some limit A and that for large r. 



Letting S 



A(r) - A 
completes the proof. 



< r 



-1/2+5 



It remains to show that the above A is equal to Aoo, the spectral gap of 
the dynamics on the infinite-volume lattice. Let {at) denote the Glauber 
dynamics on the infinite volume lattice, let (c^) be the dynamics starting 
from the all-plus configuration and define 

6=P(CT+(0) = 1)-P(CT,-(0) = 1) , 
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which in the special case of no external field is simply 2 (P ((7^(o) = l) — • 

Claim 4.4. The above defined o-nd A satisfy — exp(— A)| = 0(^^), 
and in particular lim^^oo Ct^* = exp(— A). 

Proof. Put n = exp(t) and let {X^) and (X^") denote the Glauber dynamics 
starting from the all-plus and all-minus configurations respectively on the 
torus with side-length r = 3 log^ n. Consider the monotone coupling of 
(Xf^) and (X^). By the symmetry of the torus, the expected number of 
disagreements at time t is given by r"^ (F{X^{o) = 1) — P(X^^(o) = 1)) and 
hence 

F{X+{o) = 1) - r{Xf{o) = 1) < 2max||P^o(Xt € •) " /^IItv 

XQ 

< 2r^ (P(X+(o) = 1) - P(Xr(o) = 1)) , 

(4.8) 

where the first inequality is by definition of the total variation distance. 

Next, identify the vertices of with those in {x € : ||a;||oo < '''/2} 
and couple {X^') and {at) via identical updates to identified vertices. By 
another simple application of the disagreement percolation argument (as in 



the proof of Lemma 3.6), 

P(a+(o)/X+(o)) <n-io<^. 
Combining this with ( |4.8[ ), 

r-Vax||P^„(Xt G .)-^\\^^-n-^^'^ < it < 2 max ||P^„(Xt G •)-MllTv+n-^°'^. 



Applying (4.1) and (4.3) for a choice of s = (lOd/a*) loglogn now gives 

it < 0(n^^'^) -I- 2r'^/^e^'^^''^*^*^'*-*+'^"'°^'°§" < e--'*(^)*+'=i l°glog" 

and 

C > ^^Q-Mr)t _ ^-lOd > -A(r)i-ci loglogn 

for some constant ci > 0. Altogether, as r = 31og'^ n and t = logn. 



log 



exp(-A(3i3)) 



log log n log t 
< ci ^ = ci - 



which combined with Lemma |4.3| completes the proof 



t 



A result of HoUey [1^ shows that lim^^oo ^* = exp(— Aqo)- We will 
show how this result is quickly recovered from our proof. Plugging the 
test-function /(cr) = l{o-(o)=i} — '^^{a{o)=i} iiito the characterization of the 
spectral gap as the slowest rate of exponential decay of the semigroup gives 

exp(-Aoo) > lim i]/^ = exp(-A) . 

t— >oo 
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Now fix £ > and recall the Dirichlet form (2.1), according to which 



X„o = inf 



/eL2({±i}^^/.oo) Var^^(/) ' 

where fioo is the stationary measure of the infinite-volume Ising model. For 
any / G L^({ibl}^ , /ioo) with S^^^ (/, /) we can find a sequence of functions 



fn G L^({±1}^'', /ioo) each of which depends only on a finite number of spins 
such that fn_^f in L^{{±1}'^ ,/Uoo) and £f,^{fn, fn) ^ ^i^^UJ) (see e.g. 



the proof of 24, Lemma 4.3]). So take g E L^({ibl}^'', /ioo) depending only 



on a finite number of spins such that 

Var^^(g) " Var^^(/) ' 

For some large enough M we have that g is a function of the spins in the 
box Bm = {x gZ*^ : \\x\\oo < M}. We compare the Glauber dynamics on 
Z'^ to that on for some large r, identifying the vertices of the latter with 
those in {a; G : ||x||oo < i"/"^} and denoting its stationary distribution 
by Hr- By the strong spatial mixing property, the projection of fir on Bm 
converges to the projection of /ioo on Bm as r ^ oo (in fact, an assumption 
weaker than strong spatial mixing would already infer this, e.g. uniqueness), 
hence 

—- = hm > hm A(r = A , 

Var^^(5) r^oo \ai f,^{g) r^oo 

where the inequality follows from the characterization of the spectral gap 

by the Dirichlet form. This implies that 

Aoo=inf f;°°^Y,| >A-£, 
/ Var^^(/) 

and letting e ^ gives that Aoo = A, as required. This completes the proof 
of Theorems [T] and |2j ■ 



5. Cutoff for other spin-systems 

While the proof above was given for the ferromagnetic Ising model, its 
arguments naturally extend to other well-studied spin-system models. In 
this section we outline the minor modifications one needs to make in order 
to obtain cutoff for general monotone and anti-monotone systems, thereby 
proving Theorems [5] and [6| 

A key prerequisite for our proofs is a uniformly bounded log-Sobolev con- 
stant for the dynamics. This was proved in great generality in [2^ for 
Glauber dynamics on spin-system models on the lattice where strong spa- 
tial mixing holds, and that result carries to periodic boundary with minor 
adjustments. 
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Otherwise, the only specific property of the Ising model used in our argu- 
ments is the monotonicity of the model. The Glauber dynamics (1.2) for a 
spin- system is monotone (or attractive) if there is a partial ordering of the 
state space ^ according to which the transition rates satisfy 

a{x)c{x, a) < r](x)c{x,r]) for all 77 ^ o" with r]{x) = (7{x) . (5.1) 

Under this condition, the system allows a monotone coupling of the Glauber 
dynamics. We use the monotonicity of the model in precisely two locations. 



(1) In Lemma 3.9 we prove that the update support is sparse w.h.p. by 
showing that in most of the tori that make up the barrier dynamics the 
all-plus and all-minus configurations couple (hence the final configura- 
tion is independent of the projection of the starting configuration onto 
those tori). This is a crucial application of the system's monotonicity. 

(2) In equation (4.8 ) of Claim 4.4 we relate the expected number of disagree- 
ments in the monotone coupling of the dynamics to the total-variation 
distance of the dynamics from stationarity. 

The arguments in both of these steps hold essentially unchanged for any 
monotone dynamics. 

An anti-monotone system is one where the reverse inequality in (5.1) 
always holds. When the underlying geometry is a bipartite graph, there is a 
standard transformation of an anti-monotone system into a monotone one. 
Let V = Ve L) Vo he a partition of the sites such that there are no edges 
within Ve or Vq. We define a new partial ordering :<* on {±1}^ as follows: 
For two configurations a, rj we have a ^* if c^v) < r]{v) for all f G T4 and 
a{v) > r]{v) for all v G Vo- It is easy to verify that an anti-monotone system 
under the standard partial ordering ^ is a monotone system under :<* . 

To derive Theorems [5] and |6j note that the d-dimensional lattice with 
periodic boundary conditions (Z/rZ)'^ is bipartite if and only if r is even. 
Crucially, our proof for cutoff on (Z/nZ)*^ only required the dynamics to be 
monotone on the smaller tori (Z/rZ)'^ for r = r{n) as defined above (3.1). 
As we noted there, we have freedom for our choice of r (basically any choice 
between log^"*""^ n and log*^^"^^ n would do) and in particular we can let r be 
even. Applying the above transformation therefore establishes cutoff for the 
anti-ferromagnetic Ising model and the gas hard-core model. 

As a side note, in our companion paper [26] we establish that any model 
with soft interactions (e.g. the Potts model) has cutoff at high enough tem- 
peratures. One of the key challenges there is to show that the update support 
is typically sparse. We further address non-periodic boundary conditions, 
which break down the symmetric structure of the torus, thereby letting the 
smaller tori (which we analyze in our L^-L^ reduction) play different roles 
in the dynamics according to their vicinity to the boundary. 
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